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Abstract 



A confining, Goldstone theorem preserving, separable Ansatz for the lad- 
der kernel of the two-body Bethe-Salpeter equation is constructed from phe- 
nomenologically efficacious u, d and s dressed-quark propagators. The sim- 
plicity of the approach is its merit. It provides a good description of the 
ground-state isovector-pseudoscalar, vector and axial-vector meson spectrum; 
facilitates an exploration of the relative importance of various components 
of the two-body Bethe-Salpeter amplitudes, showing that sub-leading Dirac 
components are quantitatively important in the isovector-pseudoscalar meson 
channels; and allows a scrutiny of the domain of applicability of ladder trunca- 
tion studies. A colour-antitriplet diquark spectrum is obtained. Shortcomings 
of separable Ansatze and the ladder kernel are highlighted. 

Pacs Numbers: ll.10.St, 14.40.-n, 24.85. +p, 12.40.Yx 
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I. INTRODUCTION 



The spectroscopy of light-quark mesons is made interesting because of the role played by 
dynamical chiral symmetry breaking, the natural scale of which is commensurate with other 
scales in this sector. It also explores quark and gluon confinement because most vector and 
axial-vector mesons have masses that are more than twice as large as typical constituent- 
quark masses. Covariant, constituent-quark potential models are a useful tool in the study 
of this problem |]T|. 

A salient feature of the strong interaction spectrum is the fact that m 2 — m 2 pa 30m 2 , 
which may be compared with the vector-pseudovector splitting m 2 ai — m 2 pa 1.7m 2 . Fur- 
thermore, the pion mass must vanish in the chiral limit; i.e., when the current-quark mass 
vanishes, whereas m 2 ai —>■ m 2 . (We note that a vanishing current-quark mass does not en- 
tail a vanishing of the constituent-quark mass.) These observations are an indication that 
the Goldstone-boson character of the pion is a particular and crucial feature of the strong- 
interaction spectrum. That these features are difficult to capture in potential models is well 
illustrated in Refs. |] 

An efficacious framework for studying meson spectroscopy is provided by the QCD 
Dyson-Schwinger equations [DSEs] , which include the "QCD gap-equation" (quark DSE), 
that has proven useful in the study of quark confinement and dynamical chiral symmetry 
breaking, and the covariant, two-body bound-state Bethe-Salpeter equations [BSEs]. With 
one exception, Ref. [fj, all spectroscopic studies to date have employed the rainbow- ladder 
truncation of the quark DSE and two-body BSE, which is defined as follows. Rainbow ap- 
proximation specifies that the dressed quark-gluon vertex in the quark DSE is replaced by 
the bare vertex: T a (k,p) = 7 A( A a /2, where {A a }^ =1 are the colour Gell-Mann matrices. This 
equation is then solved with a given model form of the dressed-gluon propagator, D^k), to 
yield a dressed-quark propagator, S(p) = l/[*7 ■ pA(p 2 ) + B{p 2 )]. The kernel of the ladder- 
approximation to the two-body BSE is then the customary ladder kernel but with D^ v {k) 
and S(p) employed in place of free-particle propagators. 

It has been shown [|J that for any D lll ,{k) that leads to the dynamical generation of a 
fermion mass in the chiral limit; i.e., to dynamical chiral symmetry breaking, the isovector- 
pseudoscalar meson BSE necessarily admits a P 2 = bound-state solution (P M is the total- 
momentum of the dressed-quark, antiquark system). No fine-tuning is necessary to ensure 
this outcome and one thus has a natural understanding of the pion as both a Goldstone boson 
and a bound-state of a strongly-dressed quark and antiquark. This outcome is the result of 
an equivalence, in the chiral limit, between the quark DSE and the isovector, pseudoscalar 
meson BSE. This equivalence is an intrinsic feature of the DSEs, which persists in more 
sophisticated truncation schemes [§,§J. 

The most extensive and phenomenologically successful spectroscopic studies in the 
rainbow-ladder framework are those of Ref. |J, in which the quark DSE is solved numer- 
ically for spacelike-p 2 using a model gluon propagator. In Landau gauge the behaviour of 
the gluon propagator is constrained by perturbation theory for k 2 > 1 — 2 GeV 2 [0] and one 
models the infrared behaviour, which is presently unknown. Such studies have the ability to 
unify many observables via the few parameters that characterise the behaviour of the model 
dressed-gluon propagator in the infrared. 

In this approach, solving the meson BSEs is complicated by the fact that these equations 
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sample the dressed- quark propagator off the spacelike-p 2 axis. In Ref. || this difficulty was 
circumvented by employing a derivative expansion of the dressed-quark propagator functions, 
A(p 2 ) and B(p 2 ), and estimating the error introduced thereby This, however, obscures the 
discussion and exploration of the role of quark and gluon confinement, a sufficient condition 
for which is the absence of a Lehmann representation for the dressed-quark and dressed-gluon 
propagators. The problem becomes more acute in studies of scattering observables. 

An algebraic parametrisation of a confining dressed-quark propagator, based on numer- 
ical solutions of model quark-DSEs, has been used successfully in studies of a large range 
of mesonic scattering observables; for example: f n , r n , the tc-ti scattering-lengths and the 
pion electromagnetic form factor |J; fx and the charged and neutral kaon electromagnetic 
form factors |§; the anomalous 77? — > 7 |TI| and 77? — > Tin transition form factors JTTJ. In 



these studies the dressed-gluon propagator is only specified implicitly insofar as the model 
dressed-quark propagator can be used as a constraint on its form via the quark DSE. 

It would be useful to make this connection explicit. However, given a dressed-quark 
propagator it is not possible, in principle, to invert the quark DSE and extract a dressed- 
gluon propagator; one reason being that the quark DSE involves the dressed-quark-gluon 
vertex, which depends implicitly on both the dressed-quark and dressed-gluon propagators. 
In the peculiar case of the rainbow truncation this particular difficulty, at least, is eliminated. 

It is known that the rainbow truncation is only quantitatively and qualitatively reliable 
in Landau gauge [|12|], which means that it is inappropriate to infer a connection between 
a given phenomenologically-constrained model dressed-gluon propagator and a solution of 
the gluon DSE, such as those obtained in Refs. [0,0' i n an y other covariant gauge. The 
quark DSE is, in general, a pair of coupled, non-linear integral equations for A(p 2 ) and 
B{p 2 ). In Landau gauge, the kernel in the equation for A{p 2 ) is sufficiently complicated, 
even in rainbow truncation, that it is not possible to invert the equation without introducing 
kinematic singularities. An explicit connection between the dressed-quark propagator and 
a confining dressed-gluon propagator via the inversion of the quark DSE is therefore not 
possible. 



A goal of this study, and another [14]], is to explore the extent to which pion and kaon 



observables, as embodied in the model dressed-quark propagators employed in Refs. PHTTf, 
constrain the properties of all light-quark mesons. As we have described, it is not possible 
to explicitly construct a dressed-gluon propagator from these dressed-quark propagators. 
However, one may adopt a purely phenomenological approach in order to construct a simple, 
confining Ansatz for the kernel of the two-body BSEs that is constrained by the pion and 
kaon scattering observables. This facilitates the present exploration of the extent to which 
a confining, Goldstone theorem preserving, BSE approach can generate the ground state 
spectrum of light-quark mesons. This approach allows one to easily identify those channels to 
which it is applicable and those in which it is inadequate, and to explore the extent to which 
observable properties are influenced by sub-leading Dirac components in the Bethe-Salpeter 
amplitude; e.g., the influence of the pseudovector, 757 ■ P, piece of the pion Bethe-Salpeter 
amplitude on the pion mass and decay constant. Such terms have been neglected in almost 
all studies undertaken to date. A similar situation holds for the calculation of hadronic 



coupling constants and associated form factors for processes such as p — > nn [|T^,|TH] and 



p — > 77T [17||. Semi-phenomenological qq Bethe-Salpeter amplitudes for the leading Dirac 



covariant are currently used to facilitate the necessary integrations. The amplitudes provided 
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here have a significant amount of dynamical justification and yet are simple enough to allow 
a more realistic study of hadron couplings. We also take this opportunity to explore diquark 
correlations. 

It is with this goal in mind that, in Sec. |I|, we construct a constrained, confining, flavour- 
dependent, separable Ansatz for the dressed-ladder kernel of the two-body BSEs. Light- 



quark meson spectroscopy is discussed in Sec. |ITJ. In Sec. JV] we consider the spectroscopy 
of colour-antitriplet quark-quark (diquark) correlations, which are bound in dressed-ladder 
truncation. This is a defect of the truncation, which is due to the fact that there are no 
repulsive terms in the kernel at this level of truncation. It is a peculiarity; repulsive terms 
appear at every higher order, with "order" referring to the number of explicit dressed- 
gluon propagators in the kernel, and these eliminate the diquark bound states |§. Diquark 
spectroscopy is nevertheless of contemporary interest because a number of studies of the 
nucleon Fadde'ev equation have proceeded under the assumption that the quark-quark T- 
matrix can be represented as a sum of simple diquark-pole terms, and that the mass splittings 
are such that only the lowest mass poles need be retained in solving the reduced two-body 
equation that results plfnj. We summarise and conclude in Sec. M. 



II. SEPARABLE ANSATZ FOR THE BETHE-SALPETER KERNEL 

The Dyson-Schwinger equation for the Euclidean-space dressed-quark propagator 
(Schwinger function), 

S(p) = -H -pMp 2 ) + «s(p 2 ) = n . pA{p l ) + B{f) , (i) 

can be written as 

S-\p)=^l■p + m + J0^g 2 D, u (( P -kr) 1 ^S(k)Tl(k,p) (2) 

where m is the (bare) current-quark mass. The Euclidean Dirac matrices satisfy the algebra 
{linlv} = 25^, where 5^ v is the Kronecker-delta, and a - b = J2t=i a ibi- In Eq. fl2|), D^ u (k) 
is the dressed-gluon propagator and T At (/c,p) is the dressed quark-gluon vertex. 
The homogeneous BSE for a quark- ant iquark bound state is 

T rs (k;P) = J ^K^( q , P] P)[S f A q + ^P)T( q] P)S h ( q -(l-OP)] tU , (3) 

where P is the centre-of-mass momentum; the f\ -flavour- quark carries momentum P f t = 
q + £P and the / 2 -flavour-antiquark momentum is pj 2 = — q + (1 — £)P; and K( q ,p; P) is 
the quark- ant iquark scattering kernel. 

The ladder approximation is defined by: 

Y J Iroto Iu F s F ( y J lu D s D , (4) 

/ r c t c V / u c s c 

where {F, C, D} indicate flavour, colour and Dirac indices. 
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In rainbow approximation; i.e., using 



r^,p) = 7 ,y (5) 

in Eq. (g), then, if the dressed-quark propagator is known, Eq. (0) can be used to constrain 
an Ansatz for the ladder approximation to the kernel of Eq. (^). 

Herein, solely because of its inherent simplicity, we employ a model of the form 

g 2 D^(p-k)=5^A(p-k) . (6) 

In being proportional to 5^, this has the appearance of a Feynman gauge propagator. The 
appearance is misleading, however. A fundamental Slavnov- Taylor identity in QCD entails 
that the longitudinal piece of the dressed-gluon propagator must be independent of interac- 
tions. Equation @ has the transverse and longitudinal components dressed in exactly the 
same way. A propagator of this form could only arise if the gauge parameter was chosen so 
as to completely cancel the transverse interaction contributions; i.e., if the gauge parameter 
dependent, longitudinal piece of the gluon propagator is interaction dependent. Therefore 
Eq. (|H) can only provide a model effective-potential: A(p — k), defined in this way, cannot in 
principle be related to solutions obtained in studies of the gluon DSE, such as Refs. |T|,|T3 . 
We will describe Eq. (Q) as Feynman-like gauge. 
One can write, without loss of generality, 

oo 

A(p-k) = Y J ^n(p\k 2 )p n k n ±U n (p-k) , (7) 

n=0 

where p is the unit-magnitude direction- vector for p and {U n (x)}^ =0 are the complete set of 
orthonormal Tschebyshev functions, which satisfy 

- f 1 dx Vl-x 2 Ui{x) UAx) = Sij . (8) 

7T J-l 

Translational invariance is preserved if all contributing Tschebyshev moments are retained. 

The quark DSE, Eq. (§), represents two coupled, non-linear integral equations for A(s) 
and -B(s), where s = p 2 . In rainbow approximation, Eq. (||), and using Eq. (|), these 
equations are 

At?) = f + | / A((p - *)')„■ k k2A{k ^l ] B{k2)2 , 0) 

Bi " 2) = m + f / (0 A < (p - fc)2 > m^TWr ■ m 

The simplicity inherent in Feynman-like gauge is obvious. 

Introducing the Tschebyshev expansion for A(p — k), Eq. (|7|), these equations become 

A{s) = 1 + M^Jo dt * A i(s,t)(Tv(t) , (11) 
1 r°° 

B(s)=m + — dttA (s,t)a s (t) , (12) 

OTl z JO 
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from which one observes that, in rainbow approximation and in Feynman-like gauge, the 
quark DSE is only sensitive to the zeroth and first Tschebyshev moments of A(p — k). Hence, 
translational invariance of the kernel of the quark DSE is not lost as long as the zeroth and 
first Tschebyshev moments are retained. 

A constrained kernel can now be obtained by employing a rank-N, separable approxima- 
tion for the Tschebyshev moments: 



A' 



A n (M) = £^(s)^(t) 



(13) 



i=i 



The simplest such approximation is rank-1, which is considered herein; i.e., one writes: 

F \s) = G(s) = I (B(s) - m) , F^s) = F(s) = ± (A(s) - 1) , (14) 



where a and b are fixed constants, which are to be determined, and A(s) and B(s) are the 
functions that appear in the quark propagator. As will be seen below, this particular choice 
for Fq 1 and is sufficient to ensure that Goldstone's theorem is preserved. 

Substituting Eqs. ([14]) via Eq. (|i~3|) into Eqs. flllD and (|12|) one finds that this latter pair 
of equations [i.e., the quark DSE] is solved if, and only if, 



1 r°° 

— J Q dtt 2 [A(t)-l]a v (t) , 
1 r°° 

— dtt[B(t)-m]a s (t) . 

7T JO 



(15) 
(16) 



One now has a rank-1, separable Ansatz for the kernel of the BSE for like-quarks, which 
is completely determined by the propagator of that quark; i.e., Eq. (H) with 



(17) 



g 2 D^{p -k) = 5^A(p — k) = 5, u G{p 2 ) G{k 2 ) + p -kF{p 2 ) F{k 



With u- and <i-quarks treated as indistinguishable, except for their electric charge, Eq. (Jlj 
can be used in the study of the BSE for tt, uo and p mesons, for example. 

A simple generalisation of this Ansatz to meson-like bound states with arbitrary flavour 
content is obtained via the identification 



S h (k + £P) A(p - k) S f2 (k - (1 - = 
S h (k + £P) {i [G h {p 2 ) G h {k 2 ) + G h {k 2 ) G h (p 2 
+p-k\ [F h (p 2 ) F h ik 2 ) + F h (k 2 ) F h (p 2 )} } Sf2 (k-(l- OP) 



(18) 



wherever it appears in the kernel of a given BSE. This can be used in the study of the BSE 
for K and K* mesons, for example. 

We observe that, once the propagators for quarks of flavours fi and f2 are known, Eq. ([18]) 
provides a constrained, separable Ansatz for the kernel oi the Bethe-Salpeter equation. If the 
dressed-quark propagators have no Lehmann representation then this kernel is free of quark 
and gluon production thresholds and may therefore be described as confining. As remarked 
above, this Ansatz for the kernel is not equivalent to an Ansatz for the gluon propagator and 
it is inappropriate to infer comparisons with solutions obtained in studies of the gluon DSE, 
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such as Refs. [f7yi3|. Such comparisons can only be made when one employs a gauge- fixing 



procedure that does not violate the relevant Slavnov- Taylor identity; for example, Ref. [20 



which employs Landau gauge and is not separable. We note that any attempt to construct 
a constrained, separable Ansatz in other than Feynman-like gauge will introduce kinematic 
singularities in the analogue of Eq. (p~8|) . 

The BSE is solved in the rest frame by setting P = (0, 0, 0, iM) in Eq. (|3|) and details of 
this for the case of identical quarks (fi = ^2) are given in Appendix [AJ The generalisation 
to meson-like bound states with arbitrary flavour content is straightforward using Eq. (|18|) . 
In general, the ladder truncation of the BSE reduces to a finite matrix equation that admits 
solutions for discrete values of the meson mass, M. 



A. Meson Decay Constant 

The canonical normalisation of the Bethe-Salpeter amplitude T is given by |^TJ 

2P» = N C J I tr D [T(k, -P)d?S h (k + £P)F(k, P)S f2 (k - (1 - OP)} (19) 



+tr D 



(27T) 

T(k, -P)S fl (k + tP)F(k, P)<%S h {k - (1 - OP) 



where T(k, P) T = C _1 r(— k, P)C defines the corresponding anti-meson amplitude. In ladder 
approximation the kernel of the Bethe-Salpeter equation is independent of the centre-of-mass 
momentum, P, hence there is no contribution of the type dK/dP to the normalisation. 
The pseudoscalar meson decay constant, f-p, is defined by: 

_ A r 

<0|tf (0) 7M 75-2-*(0)|$(P)> = P,fv, (20) 

where |$(P)) is the pseudoscalar meson state vector, A v are matrices acting in flavour 
space and \1/ is a colour triplet and flavour multiplet of Dirac spinors. For the K~ meson, 
for example, the relevant flavour matrix is, with {Aj}f =1 the Gell-Mann matrices, 



A X ~ = ^(A4 + ^A 5 ) 



/0 v/2 

I , (21) 
\0 



which gives (0|^ U (0)7 M 7 5 * S (0)|<I> X -(P)) = ^JIP^Jk- . Thus, the decay constants for the 
pseudoscalar meson solutions to the BSE given in Eq. (|3|) are defined by 

V2P,f M = (01^(0)7^75^ (0)|$a/(P)> • (22) 

To obtain an expression in terms of the Bethe-Salpeter amplitude, we note that the 
unamputated BS wave-function 

X (p, P) = S fl (p + £P)r(p, P)S h (p - (1 - op), (23) 
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can be expressed as 

rf 4 xrf 4 7/e- iP[( ^ +(1 - ?)s/] e- i( ^- p ^(0|^ /l (x)^ /2 (y)|$(P)). (24) 



(For a colour singlet bound state, xip'jP) is diagonal in colour-space.) Multiplying both 
sides by 757 • P, taking the matrix trace throughout, evaluating the integrals over p and q 
and using Eq. (|22|) one obtains 

N r d A n 

P 2 fM = J ^tr D [757 " PS fl (p + £P)I>, P)S h (P + (1 - OP)} , (25) 

which provides the relation between the Bethe-Salpeter amplitude and the canonically de- 
fined meson decay constant. In this equation V is normalised according to Eq. fliTf). 



B. Dressed Quark Propagators. 

The separable Ansatz is completely defined once the quark propagators are specified. 
Following Ref . || , the scalar and vector parts of the quark propagators are defined in terms 
of dimensionless functions: 

af/{s) = 5d*v {x) ' = 7m** 3 {x) ' (26) 



with s = p 2 , x = s/(2D), D is a mass-scale parameter, and where (A = 10 4 ): 

4(x) = — —^—2 (l - e-^iA + 
bK ' x + mj V ) 

1 _ e ~b(x 1 _ -b{x / 1 _ e -Ax \ 

— x 4^\ ba + 4 ^-)' (27) 

and 



2(x + m})-l + e~ 2ix+m f ) 
2(x + mj) 2 



Here m/ = m//v2~D- In this work the u and d quarks are considered to be identical, except 
for their electric charge. 

The dressed-quark propagator described by Eqs. ( pTf ) and ( p8|) is an entire function in the 
finite complex p 2 -plane and may therefore be interpreted as describing a confined particle . 
The ~ e _x form that ensures this is suggested by the algebraic solution of the model DSE 



studied in Ref. [p2|1 , which employed a confining model dressed-gluon propagator and dressed 
quark-gluon vertex. Furthermore, the behaviour of Eqs. (pTD and fl28| ) on the spacelike-p 2 



axis is such that, neglecting ln[p 2 ] corrections associated with the anomalous dimension of the 
dressed-quark propagator in QCD, which are quantitatively unimportant herein, asymptotic 
freedom is manifest. In Eq. (p7|) the term ~ 1/x 2 allows for the representation of dynamical 
chiral symmetry breaking and the ~ m/x term represents explicit chiral symmetry breaking. 



S 



In Ref. H the five parameters {rh u , 6q, . . . , 63} in Eqs. (|27j ) and fl28|) were varied in 
order to determine whether this model form could provide a good description of the pion 
observables: /„-; m n ; (qq); r n ; the n-n scattering lengths and partial wave amplitudes; and 
the electromagnetic pion form factor. A very good fit was found with the w-quark parameter 
values listed in Eq. (p^): 

u- quark s- quark 

(29) 



771/ 


0.00897 


0.224 


bl 


0.131 


0.105 


bl 


2.90 


2.90 


b{ 


0.603 


0.740 


bl 


0.185 


0.185 



The scale is set with D = 0.160 GeV . This same model also provides a good description of 



the 7*7r — > 7 |T0[ and "fir* — > ttti |TT[ transition form factors 



Dyson-Schwinger equation studies |2I] indicate that while it is a good approximation 
to represent the u- and <i-quarks by the same propagator, this is not true for the s-quark. 
For example; contemporary theoretical studies suggest that 2m s /{m u + m^) ~ 17 — 25 f25| 
and (ss) ~ 0.5 — 0.8 (uu) ||26|| , which is a nonperturbative difference. In Ref. ||, with this 
in mind, the model forms in Eqs. ( |2"T| ) and fl2"E| ) were employed in a study of the kaon 
observables: fx] (ss); r K o; r K ±; and the electromagnetic form factors of the charged and 
neutral kaon. The sensitivity of these observables to rh s and (ss) was too weak for an 
independent determination and therefore fh s = 25rh u and 6q = 0.869; which ensures (ss) = 
0.8(uu), were chosen for consistency with other theoretical estimates. The parameter 
was allowed to vary to provide a minimal residual difference between the u/d- and s-quark 
propagators and a very good fit to the kaon observables was obtained with the value listed 
in Eq. Q2|). 

To complete the specification of the constrained separable approximation to the kernel of 
the Bethe-Salpeter equation, the quantities a and b in Eqs. flT5| ) and flTBj ) must be determined. 



However, using Eqs. (^7]) and (|28| ) neither a nor b is finite. Equations ( |I~5"D and (|1^) only 



yield finite values if the large spacelike-x behaviour of ay and a$ is such that: 

Mx) = + O (-^t) , <ts(x) = + O (-^t) , (30) 

x + m 2 \x 2+s J x + m 2 \x 2+s J 

for any 5 > 0. Dynamical chiral symmetry breaking in QCD entails that at large x [up to 
corrections ~ [lnx] -7 , 7 < 1] as( x ) — fn/(x + rh 2 ) + 0(x~ 2 ) and hence no quark propagator 
that properly incorporates the momentum dependence at large-x due to dynamical chiral 
symmetry breaking will yield finite values of a and b. (This behaviour is tied to the necessary 
divergence of the quark condensates in QCD; necessary because condensates are related to 
two-point Schwinger functions evaluated at zero relative Euclidean spatial separation.) 

To complete the specification of the constrained, separable Ansatz one must therefore 
incorporate an ultraviolet regularisation in the propagator: 

6 w x + m 2 f V / 

1 - e- b ^ x 1 - e~ b i x ( j f 1 - e - biX \ 1 - e -( £ s x ) 2 

bo+b 2 7 f (31) 



b[x b{x \ b 4 x 



6cX 
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which introduces three new parameters: ey, eg, e|, that are not determined by the studies 
of Ref. 0. The parameter ey = 0.1 is chosen so as to ensure that a v Keg are numerically 
good approximations to a v on the domain < x < 3; our results are not sensitive to the 
domain x > 3. It is not varied but we have established that our results are insensitive to 
it; i.e., that changes can be absorbed into a change in e$. The regularisation parameters 
modify the large-p 2 behaviour of the propagator, which entails that the light-quark mass 
values must be re-fit {m q — > m q ). 

Equations (|T^-|TB|), (P"8|), (|2DJ), fl3"T|) and fl3"2| ) completely specify the constrained, confining, 
separable Ansatz for the ladder-kernel of the Bethe-Salpeter equation. The numerical studies 
proceed by varying the four parameters m, and e f s in order to fit U /K and m n/K and then 
predicting the ground state spectrum of octet mesons. Diquark systems are also studied. 

III. MESONS 

The Bethe-Salpeter equation considered for a bound state of a quark of flavour fi and 
an antiquark of flavour f 2 is 

r(p, P) = -\f j^y A (p - ?bA (<? + £P)r( q , P)s f2 ( q - (1 - £)p) v (33) 

In this equation A(p — q) is obtained from Eq. flTSP with Sf t obtained from Eqs. ([H]) and 
(|32|) using the parameters in Eq. ( p9| ) and Table |. This equation is solved in each channel 
as an eigenvalue problem of the form KT = X(P 2 )T, with the bound state mass identified 
from A(P 2 = — M 2 ) = 1. 

A. Scalar and Pseudoscalar Mesons. 

1. h = u/d = h 

In this case the requirement of charge conjugation invariance for the neutral mesons 
entails £ = 1/2. The form of the charge parity C = ±, fi = u/d = f 2 Bethe-Salpeter ampli- 
tudes, Tq, obtained as solutions at the mass-shell point P 2 = —M 2 are given in Eqs. ( |B1| - ]B4| ). 

These equations expose a shortcoming of separable Ansatze: the pseudoscalar and pseu- 
dovector pieces of the pseudoscalar Bethe-Salpeter amplitude are characterised by the same 
function, which is not the case in general. 

The calculated eigen- vectors are given in Eq. ( |B6| ) and the bound state masses in Table ||. 

The separable Ansatz for the kernel of the Bethe-Salpeter equation yields mo-+ = 
when fhu/d = 0. This is a necessary consequence of the equivalence between the isovector- 
pseudoscalar BSE and the quark DSE in this chiral limit ||, which is preserved in the 
approach described herein and discussed in detail in Ref. [[J. 

One might be tempted to conclude from Eq. (|B6|) that for C = + states the leading 
Dirac component of the amplitude dominates; i.e., the pure 75 component dominates for 
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the _+ state and the pure Id component for the ++ state. Indeed, it is an often used 
approximation to neglect sub-leading Dirac components of the Bethe-Salpeter amplitude 
in ground state studies using the Bethe-Salpeter equation. Considering Tables [TT] and fT| 
one observes that while this is a good approximation for the heavy ++ state, it represents 
an erroneous conclusion for the light CT + state, for which the sub-leading, axial-vector 
component provides 17% of the mass and 39% of the decay constant. This feature is also 
seen in Ref. M. 

Table |I| shows that the separable Ansatz for the BSE kernel yields a large ttiq++ — mo-+ 
splitting without fine tuning, thus reproducing this characteristic feature of the strong inter- 
action spectrum. The ++ state can be identified with the ao(980) meson. The discrepancy 
between the calculated and observed masses is consistent with the contention that this state 
involves a considerable K-K admixture, which can be represented as a contribution to the 
Bethe-Salpeter kernel but is absent in ladder approximation. 

No J PC = + ~ states have been observed in the strong interaction spectrum. However, in 
general, as observed in Ref. p^fl , the Bethe-Salpeter equation admits solutions of this type. 



The amplitudes for such solutions characteristically differ from their C = + counterparts 
by the factor p ■ P which is odd under charge conjugation. Such states have no analogue 
in quantum mechanics since, for equal-mass constituent-particles on shell, p • P — 0. In 
this context one observes that our separable Ansatz for the kernel of the Bethe-Salpeter 
equation yields very heavy and + ~ states, with m — ~ 10m -+ and m +- ~ 2m ++; 
i.e., it yields results consistent with the observed strong interaction spectrum. 



2. h = u/d, f 2 = s 

The Bethe-Salpeter equation for u-s states is Eq. ( |33| ) with f\ = u/d, f 2 = s. Consider 
first the pseudoscalar (kaon) channel. A value of £ is determined by ensuring that the electric 
charge of the K° is zero in impulse approximation [[|. The value obtained in Ref. || with 
empirical Bethe-Salpeter amplitudes is £ = 0.49 (~ 0.5), while the present work requires 
£ = 0.56 (~ 0.5). There is only a weak sensitivity of masses to changes in £ of this magnitude 
throughout this work. 

The solution amplitude is given in Eq. (|B7| ) with the calculated eigen-vector given in 
Eq. ( |B8| ) and the mass in Table |T[ We note that there is no charge-parity, C, symmetry for 
bound states of distinguishable quarks. 

The leading, pseudoscalar Dirac amplitudes (Xi u , Ai s ) again appear to be dominant 
for the kaon, however, as for the pion, the sub-leading, axial-vector amplitudes (A 3u , A 3s ) 
contribute significantly to the mass (17%) and decay constant (33%). 

Each type of covariant in the kaon solution is weighted by two amplitudes that de- 
scribe the internal momentum dependence in terms of functions that relate to the dressed- 
propagators of the u/d- and s-quarks. These are found to have approximately equal in- 
fluence in the solution. For example, from Eq. ( p8|) and Table | one calculates that 
^i/b u = 9.4 GeV -1 and A2/&S = 9.2 GeV -1 . This means that, using the constrained, 
separable Ansatz, the kaon Bethe-Salpeter amplitude for the pseudoscalar covariant is an 
approximately even mixture of the u- and s-quark mass functions. 

It is clear from Table | that 
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which is essentially the same as the ratio obtained from the mass values in Eq. fl29|) and is 
in the range (17-25) suggested by other theoretical analyses []2~5| . With D = 0.160 GeV 2 , 
™u/d = 0.00811 corresponds to m u / d = 4.6 MeV and rh s = 0.198 corresponds to m s = 

112 MeV. These values should not be compared directly with values of va u j~^ s quoted 
by other authors because the regularisation of the vacuum condensates employed herein, via 
the parameters ey and e£ in Eqs. ([H]) and (|32"D , is unconventional and enters through the 
quantities a and b in Eqs. (p~5f ) and (ff6|). The ratio, Eq. (^4|), is likely to be less sensitive to 
this difference and therefore provides a meaningful point of comparison. 

The dressed-quark propagators we employ are confining, with the dressed-quark mass 
being a function of p 2 , M(p 2 ), such that there is no dressed-quark mass-pole; i.e., no solution 
of the equation p 2 [A^(p 2 )} 2 + [B^(p 2 )] 2 = 0. A simple estimate of the value of the mass 
function that is most important in calculations of meson observables is obtained from the 
solution of — p 2 [Af(p 2 )} 2 + [B^ (p 2 )} 2 = 0, which might be called the Euclidean constituent- 
quark mass, M^. With the parameter values used herein, MJ. = 315 MeV and M S E = 
397 MeV. 

Our calculations yield no true + eigenstate with a mass less than 2 GeV, which we 
consider to be the upper limit for the present approach. The condition for an eigenstate 
was closest to being satisfied at a mass of 1.18 GeV. This is again consistent with a large 
mo+ —mo- splitting without fine tuning. This + state might be identified with the _ffo(1430). 
Such an identification would suggest that this state, like the a (980), has a sizeable coupling 
to other channels, which contribute to its mass; i.e., that the ladder kernel is inadequate to 
properly describe this channel. 



3. n Meson. 

Ladder approximation is inadequate to properly study the 77-7/ complex. A minimal 
extension that can dynamically couple the flavour octet and singlet channels is the inclusion 
of timelike-gluon exchange diagrams. This is not considered here. 

Instead we study 

4 r d 4 a 

r„(p,P)= 1 1 



3 J (2tt) 4 

j(cosfl P - y/2 sm6 P ) 2 A u (p - q) lfl S u (q + |P)r„(g, P)S u (q - ±P) 7m + 

§0/2 costfp + sin 6 P ) 2 A s (p - q) llx S s {q + ±P)T v (q, P)S s (q - |P) 7 J , (35) 



with 



A f (p-q)=G f (p 2 )G f (q 2 ) + P -qF f (p 2 )F f (q 2 ). (36) 
Eq. ( P^D is the projected Bethe-Salpeter equation for the meson whose flavour structure is 

F v = A 8 cos 6 P - A sin 6 P , (37) 
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with A = y2/3diag(l, 1, 1) and Op an octet-singlet mixing angle. The exact kernel of the 
Bethe-Salpeter equation would lead to a prediction for Op. 

With the kernel considered herein, Op is treated as an external parameter on which the 
mass and other properties of the ?7-meson depend. For example, in this case the expressions 
for the normalisation of the Bethe-Salpeter amplitude, Eq. (|I9~D, and the decay constant, 
Eq. (|25|), are #p-dependent. The modified forms are given in Eqs. (|B^ ) and ( P10|) , respec- 
tively. 

The form of the positive charge parity solution of Eq. (|35p is given in Eq. ( |B11| ). The 
calculated mass, decay constant and eigen- vector, at a number of values of 8p, are given in 
Eq. ( pi2|) . The experimental values of the mass and decay constant are given in Table ||. 

In this case the sub-leading Dirac amplitudes contribute ~ 14% to the mass and ~ 26% 
to the decay constant. 

The constrained separable Ansatz favours a small positive value for the mixing angle, 
Op. This can be compared with 9p = —10° estimated in Ref. [2S|. As remarked therein, 
however, there are large uncertainties in this value. 

The ?/-meson can be studied via the projection of the Bethe-Salpeter equation orthogonal 
to that in Eq. ([35]), which is obtained from this equation under 0p^0p — n/2. As remarked 
above, one expects timelike gluon exchange, forbidden in the flavour octet channel, to be 
important in this mainly singlet channel. The results in Eq. ( [B 1 2|) , which one might compare 



with the experimental values of M^ xpt ' = 958 MeV and /„/ = 89.1 ± 5 or 77.8 ± 5 MeV, 
may be interpreted as a guide to the importance of such contributions in this channel and 
emphasise the necessity to go beyond ladder approximation for the rj' state. 

B. Vector and Axial- vector Mesons 

The ladder approximation to the Bethe-Salpeter equation for vector and axial vector 
mesons is 

r„(p, P) = -\j (0i a (p - (? + tfWg, P)s h (g - (i - e)P) 7 „ ; (38) 

which is identical to Eq. (^) except that the Bethe-Salpeter amplitude carries a Lorentz 
index. On-shell vector and axial-vector bound states are transverse: 

P u T u (p,P) = 0, (39) 

which constrains the general form of the Bethe-Salpeter amplitude. 

1. h = u/d = f 2 

The most general form of the vector meson Bethe-Salpeter amplitude when using the 



separable Ansatz is given in Eq. ( B14|) and that for the axial-vector meson is given in 
Eq. (JETTED. 

These equations expose another shortcoming of separable Ansatze: the vector and axial- 
vector meson Bethe-Salpeter amplitudes are characterised by the same functions as the 
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pseudoscalar mesons, which is not true in general. In Ref. ||, for example, the vector meson 
amplitudes were found to be much narrower in momentum space. 

The ladder approximation does not distinguish between 1 = and 1 = 1, hence the 
vector channel corresponds to both the uo and p mesons. Similarly, the axial-vector channel 
corresponds to the f\ and a\ mesons. 

The calculated mass for these states is presented in Table || and the eigen- vectors in 
Eq. ( |B17| ). With pion and kaon physics used to fix the parameters of the quark propagators 
as described in Sec. |I[ these results are predictions. 

The sub-leading Dirac amplitudes contribute very little to the J = 1 meson masses. 

The relevant experimental value to compare the vector meson with is M^ xpt = 782 MeV, 
since it is known that pion loop dressing will lower the p-meson mass while having little 
effect on the cu-meson A recent study of this effect |T1J yields M w — M p = 21.0 MeV. 

These results in the u — d sector indicate that the u/d quark propagator parameters, 
previously set by pion physics, have produced a separable BSE kernel that captures the 
dominant physics for the ground state vector and axial vector channels. 



2. h = u/d, f 2 = s 

The general form of the Bethe-Salpeter amplitude for the u-s meson, which corresponds 
to the J p = 1" K*+ -meson, is given in Eq. ( |B18|) . We choose £ so as to ensure the neutrality 
of the K*° meson, which produces £ = 0.49 pa 0.5. 

The form of the amplitude for J p = 1 + , which is a nearly equal mix of i^i(1270) and 
i^i(1400), is simply 75 times this. The corresponding choice for £ yields £ = 0.50. The 
calculated masses are listed in Table H an d the eigen- vectors in Eq. (|B19| ) . 

The sub-leading Dirac amplitudes contribute little to the masses. 



3. h = s = f 2 

The J PC = 1 ss state, identified with the 0-meson, is computed in exactly the same 
manner as the u/p-meson except for the replacement of the w-quark propagator with that 
for the s-quark. The Bethe-Salpeter amplitude for the <fi and for the 1 ++ state, identified 
with the /i(1510) meson, are described in Appendix |B 2 c . 

The calculated masses are listed in Table || and the eigen- vectors in Eq. (|B20|) . The 
sub- leading Dirac amplitudes are again unimportant. 



4- J = 1 Summary 

We observe that the J = 1 meson spectrum is satisfactorily reproduced. These higher- 
mass states explore a larger domain in the complex quark-momentum plane than do the pion 
and kaon, which are used to constrain the separable Ansatz for the ladder kernel. This is 
an indication that a successful description of a subset of hadronic observables can translate 
into a uniformly good description of a broad range of phenomena, which is a feature that 
underlies many applications of this framework and emphasises the utility of studies such as 
that of Ref. 0. 



14 



IV. DIQUARK CORRELATIONS 



The derivation of the homogeneous Bethe-Salpeter equation from the inhomogeneous 
equation for the two-body T-matrix proceeds under the assumption that there exists a 
bound-state pole in the channel under consideration. In QCD, one expects that confinement 
ensures the absence of such poles in the quark-quark T-matrix and hence that there are no 
solutions to the homogeneous Bethe-Salpeter equation in any colour-antitriplet quark-quark 
(diquark) channel. This is supported by the studies of Ref. 0, which indicate, however, 
that one must proceed beyond ladder approximation to obtain this result. In ladder approx- 
imation one finds bound-state, diquark solutions. This is a defect of the truncation. 

Studies of the nucleon as a bound-state of three dressed-quarks using the covariant 



Fadde'ev equation have been undertaken [ 18| . The appearance of the pole in the ladder 



approximation to the homogeneous, quark-quark Bethe-Salpeter equation was used therein 
to simplify the three-body problem; i.e., to re-express it as an effective two-body, quark- 
diquark problem. This technique can also be said to underly the study of Ref. |19|] . Presently, 
the only justification for this Ansatz is the simplicity it introduces into the problem. 

Accepting this approach for the present it is then important to identify those diquark 
correlations that contribute significantly to a given three-body bound-state. As a guide one 
might assume that those diquarks whose mass is greater than that of the three-body bound- 
state under consideration would contribute little to the three-body ground-state mass. Such 
studies of the "u/d-diquark spectrum" have been reported in Refs. [p7j^9| . Herein we extend 
these studies to SUf(3). 

The ladder approximation to the homogeneous Bethe-Salpeter equation for a diquark 
correlation involving quarks of flavour fx and /2 is 



r g (p,P)= (40) 

where T denotes matrix transpose. The study of such correlations is simplified if one defines 

r^(p,p) = r- 3 { P ,p)c (4i) 



where C = 7274 is the charge conjugation matrix. It follows from Eq. (|40|) that this auxiliary 
amplitude satisfies 

If (P, P) = -f / -0jI A (P - *)7, S fl (q + £P) If (g, P) S h (q - (1 - £)P) 7 , . (42) 



It is immediately obvious that Eq. ([42] ) is identical to Eq. (|33[ ) but for a reduction in the 
(purely-attractive) coupling strength: 4/3 — ► 2/3. This observation in Ref. [2^] entailed the 



result that the mass of the scalar- (u — d) diquark is greater than rn^; and that of the vector 
(u — u), (u — d) and (d — d) correlations is greater than the mass of the ai(1280)-meson. 
(This result is true in an arbitrary covariant gauge and independent of the form of the gluon 
propagator. However, it is peculiar to ladder approximation. As discussed in Ref. |3[], any 
other truncation of the kernel of the Bethe-Salpeter equation introduces repulsive terms that 
eliminate the diquark pole.) 
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A. Scalar and Pseudoscalar Diquarks 



1. h = u/d = h 

To obtain the J p = + diquark solution of Eq. ([5]) one searches for the 0~ auxiliary 
amplitude solution of Eq. ([4^). The latter can be written in the form 

If (p, P) = G u (p 2 ) [\* -i\%yP]iy s , (43) 

which is identical in form to r+ scud in Eq. (|BT|). The 0~ pseudoscalar diquark solution is 
described by an auxiliary amplitude identical in form to r^ alar of Eq. ( |B2|) . The calculated 
masses are listed in Table [TV] and the eigen- vectors are given in Eq. ( |B21| ). 

One observes that the sub-leading Dirac amplitude contributes 11% to the + diquark 
mass. 

This result suggests that the + diquark-pole will provide a contribution to the truncated 
quark-quark T-matrix that is important in the type of Fadde'ev equation studies of the 
nucleon described above. The much larger mass found to be associated with the 0~ diquark 
correlation suggests that it may be neglected in such studies. 



2. fx = u/d, h = s 

The homogeneous Bethe-Salpeter equation in the u/d-s quark-quark channel can be 
written in the form: 

rjffo P) = -\j ^-A^ ~ + e^r*? (?, P)S.{q - (i - OPh„, (44) 

where the momentum partitioning parameter is £ = 0.56 ~ 0.5, as for the kaon. 

The solution of this equation that corresponds to the + diquark is identical in form 
to r pseud in Eq. (|B7|). The calculated mass is listed in Table [IV] and the eigen- vector in 
Eq. ( |B23|) . The sub-leading Dirac amplitudes contribute 11% to the + diquark mass. 
The magnitude of its mass is such that this correlation may be important in the Fadde'ev 
equation studies of the strange octet-baryons. No 0~ solution with a mass less than 2 GeV 
was found. This is in accord with our finding that in the + meson channel, there was 
insufficient attraction for a clear bound state. 

One observes that the diquark mass splitting M us — M U( i = 145 MeV. This may be 
compared with mg — m p w 250 MeV. One might infer from this that Fadde'ev equation 
studies, such as the ones described above, may yield the correct ordering and level-separation 
of the octet baryons. 



B. Vector and Axial-vector Diquarks 

1. h = u/d = h 

In ladder approximation the homogeneous Bethe-Salpeter equation for vector and axial- 
vector colour-antitriplet diquark correlations has the same form as Eq. (|40D except that the 
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Bethe-Salpeter amplitude carries a Lorentz index. It can be recast into the form of Eq. 
in the same manner. 

For the axial-vector (1 + ) diquark channel, the auxiliary amplitude T^(p,P) is identical 
in form to the vector meson amplitude in Eq. (|B14|) . The calculated mass is listed in Table [IV 



and the eigen- vector in Eq. (|B22| ). The axial- vector diquark mass is larger than that of the 
vector meson, in agreement with the argument of Ref . |29[ . However, it is comparable to the 
predicted scalar diquark mass. Hence the 1 + diquark-pole is likely to provide a contribution 
to the truncated quark-quark T-matrix, employed in the type of simplified nucleon Fadde'ev 
equation studies described above, that is comparable to that of the scalar diquark. 

In the 1~ channel, the auxiliary amplitude T 3C (p, P) has the axial-vector meson form 



given in Eq. ( B16|) . The calculated vector diquark mass is listed in Table [IV| and the eigen- 



vector in Eq. ( |B22| ). It is too massive to be of importance. 



Sub-leading Dirac amplitudes contribute little in these channels. 



2. /i = u/d, f 2 = s 

The auxiliary amplitude for the 1 + diquark has the same form as the vector meson 
amplitude in Eq. (BT3), while that for the 1~ diquark is simply 75 times this. The calculated 
masses are given in Table |V| and the eigen- vectors in Eq. ( |B24| ). 

These results suggest that the axial- vector diquark can be important in Fadde'ev equation 
studies of strange baryons, whereas the vector diquark can be neglected. Again, sub-leading 
Dirac amplitudes contribute little in these channels. 



3. h=s = f 2 

The auxiliary amplitude for the 1 + diquark is identical in form to that for the (f> meson 
while that for the 1~ diquark has the form of the axial counterpart (/i(1510)), both of 
which are described in Appendix [B 2 4 The calculated masses are given in Table [TV] and the 



eigen- vectors in Eq. (|B25|) . 

The low mass of the axial-vector diquark suggests that it can be important in Fadde'ev 
equation studies of all strangeness carrying baryons, whereas again the vector diquark can 
be neglected. The sub-leading Dirac amplitudes are unimportant in these channels. 



V. SUMMARY AND CONCLUSIONS 

We have constructed a crude, confining, separable Ansatz for the ladder kernel of the 
two-body Bethe-Salpeter equation [BSE] from the phenomenologically efficacious u/d and 
s dressed-quark propagators of Ref. ||. We have emphasised that no connection can be 
made between this crude kernel and the solution of the Dyson-Schwinger equation for the 
dressed- gluon propagator. 

A very good description of the ground-state, 577/(3), isovector-pseudoscalar, vector and 
axial-vector meson spectrum was obtained. Scalar mesons and the pseudoscalar rj — rf 
complex were poorly described and this is not unexpected in ladder approximation. Given 
the crudity of our construction, our results can be interpreted as demonstrating the reliability 
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and extent of applicability of the rainbow-ladder approximation to the quark-DSE/meson- 
BSE complex. 

We found that in the isovector-pseudoscalar meson channel the sub-leading Dirac com- 
ponents of the Bethe-Salpeter amplitude; i.e., those terms whose Dirac matrix structure is 
more than just 75, provide quantitatively important contributions to the mass (~ 15% ef- 
fects) and weak decay constant (~ 35% effects). These terms are unimportant in the vector 
and axial-vector meson channels. 

We saw that separable Ansatze have a number of shortcomings. In the pseudoscalar 
channel one finds that the 75 and 757 ■ P components of the meson Bethe-Salpeter amplitude 
are characterised by the same function, B(p 2 ), which is not true in general. One also finds 
that the dominant components in the Bethe-Salpeter amplitudes of the vector and axial- 
vector mesons are characterised by the same functions that characterise these components 
of the pseudoscalar mesons, B{j> 2 ). More sophisticated studies indicate that the vector 
meson amplitudes are narrower in momentum space. This indicates that the amplitudes we 
have obtained should be used with caution in the calculation of, for example, meson-meson 
scattering processes. 

The shortcomings notwithstanding, there are areas of study in hadronic physics for which 
the presently provided Bethe-Salpeter amplitudes have significantly greater dynamical jus- 
tification than that of currently used approximations. For example, hadronic coupling con- 
stants such as g p7rn |16| and g 17Tp [|17] have been reproduced from the qq structure of the 



mesons in terms of a single dominant Dirac covariant if the amplitude is allowed some 
phenomenological freedom. A more realistic treatment is facilitated by the present work. 

The ladder kernel also has the defect that it is purely attractive in both the colour- 
singlet q-q and colour-antitriplet q-q channels. This entails that it yields bound colour- 
antitriplet diquarks. This is a peculiarity of ladder approximation. Measuring "order" by 
the number of dressed-gluon lines in the Bethe-Salpeter kernel, ladder approximation is the 
lowest order kernel. Repulsive terms appear at every higher order. It has been shown |3| 
that in the isovector-pseudoscalar and vector meson channels, these repulsive terms are 
cancelled by attractive terms of the same order. This explains why ladder approximation is 
phenomenologically successful in these channels. In the colour-antitriplet diquark channel 
the algebra of SU C (3) entails that the repulsive terms are stronger; they are not completely 
cancelled and eliminate the diquark bound states ||. 

The artificial diquark spectrum we obtain is nevertheless of contemporary interest be- 
cause there have been a number of studies of the covariant, three-body Fadde'ev equation 
that use the existence of diquark poles in the quark-quark T-matrix to reduce this prob- 
lem to a two-body, quark-diquark bound-state problem. Our constrained, separable Ansatz 
indicates that such studies of the baryon should include SUf(3) scalar and pseudovector 
diquarks, since these are low in mass, but can neglect pseudoscalar and vector diquarks. 

The diquark results might be used in the following way. The study of Ref. || suggests 
that, even though colour-antitriplet states are not bound, one may associate an inverse 
correlation-length, M, with each channel; the "bound-state mass" providing an estimate of 
this. One might then construct a "pseudo-pole" representation of the quark-quark T-matrix 
(for example: ~ J2n a n [1 — exp(— [P 2 + M 2 })} / [P 2 + M 2 }) , which would not entail asymptotic 
(unconfined) diquark states but would provide for a simplification of the covariant, three- 
body Fadde'ev equation. 
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Finally, this study shows that in order to directly connect hadron phenomena with the 



dressed-gluon propagator, D fiu (k), one must start with a form of D fiu (k), as in Ref. |20 
Other approaches, while they may provide a useful phenomenology, efficacious in that it 
correlates many observables via few parameters, can only loosely constrain D liv {k) and 
hence the nature of the quark-quark interaction in the infrared. 
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APPENDIX A: BETHE-SALPETER EQUATION FOR EQUAL MASS QUARKS 



Here we present details of the solution of the ladder approximation to the Bethe-Salpeter 
equation in the case of equal mass quarks (j\ = j'2) using a separable Ansatz for the kernel. 
We begin with Eq. (gjj subject to Eqs. (|) and (|) with f ± = f 2 and £ = 1/2: 



r(p, P) = -\J tSi a (p - fffrA? + hPW(q, P)S(g - \p)i, • 



(27T) 



The general form of the scalar and pseudoscalar meson amplitudes is [^S 



-\ scalar 



(q,P) 



and 



gr(q 2 , P 2 , q ■ P)I + [g P (q 2 , P 2 , q ■ P)P^ + g u (q 2 , P 2 , q ■ P)u,{q) 
r pseud (g,P) = 

g 5 (q 2 , P 2 , q ■ P) l5 + \g P5 (q 2 , P 2 , q ■ P)P, + g u5 (q 2 , P 2 , q ■ P)u,{q) 



Hi*, 



where 



u>n(p, P) 



- PP^ 

p 2 + (p- P) 2 



u(p) \ = Ju(p,P) 2 , u^p) 



Uf,(p, P) 

\u(p)\ 



(Al) 



(A2) 



(A3) 



(A4) 



with P M [P 2 = —1] the direction-vector associated with P^. 

In Feynman-like gauge it follows from the Fierz identity that there is no piece proportional 
to [7 • P, 7- q]. For mesons which are even (odd) under charge conjugation, gi, g u , g 5 and 
#P5 are even (odd) functions and gp and g u s odd (even) functions of q ■ P. 
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Defining k^ = q^ + \P^ and = - one has 

P • u(q) = , k ■ u(q) = I ■ u(q) — 1 , k • u{p) = I ■ u{p) = q ■ u{p) . (A5) 



Multiplying Eq. (|A2| ) by I, 7 • P or 7 • u(p) and taking traces one projects out a set of 
coupled integral equations for the scalar meson amplitudes. Defining 

f(p,P) = gi(p,P) , W(p,P) = iMg P (p,P) , U(p,P) = \u(p)\g u (p,P) , (A6) 

where i M = \fP* , these equations take the following simple form: 
16 f d A q 



Sip) 

W(p) 
U(p) 



(2vr) 4 
d A q 



A(p - q) [T ff f(q) + T fW W(q) + T fU U(q)] 



d 4 q 



A(p - q) [T Wf f(q) + T ww W(q) + T wu U(q)} , (A7) 
A(p - q) u(p, P) ■ u(q, P) [T Uf f(q) + T uw W(q) + T uv U{q)} , 



where 



k ■ I \av\ 2 — \crs\ 



ry-iscalar 

=(k-l + 2k-Pl-P) \a v \ 2 + \a s \ 



T, 



scalar 
UU 



k-l 



u(q,Py 



u I 2 1 \* 1 2 
\&v\ + \&s\ 



(A8) 



and 



T^scalar 
1 }W 

T^scalar 
J-fU 



T^scalar 
1 WU 



Iff r = M 3(ct^ct 5 ) + 2iq ■ P$t(a* v a s ) 



■uf 



rp scalar 
1 UW 



\u(q)\ 
|u(g)| 



(A9) 



with 



CTy (fc 2 ) , CTy = CTy (Z 2 ) 



(A10) 



and similarly for (75. In these equations 



S{v* v vs) = Yi ( a v a s - cr v a* s ) . 
In terms of the functions 
f(p,P) = ig s (p,P) , W(p,P)=tMg P5 (p,P) , U(p,P) = \u(p)\g u5 (p,P) 



(All) 



(A12) 
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the equations for the pseudoscalar states have the form in Eq. (|A7|) but with the Ts replaced 

by 



and 





1 !f ~ K ' 

rppseud _ fi 

1 ww — \^ 




rppseud _ I , 

1 uu — 1 K 

V 


ry-ipseud 

fW 

ry-ipscud 


mpSCud 

rppseud 


ry-ipseud 

1 wu 


ry-ipseud 



\a v \ - \a s 

|2 i |2 



u(q,P)' 



(tv\-\<ts\, (A13) 



-M KKa 5 ) + 2 • P %a* v a s ) , 

|u(g)| (A14) 

n .q-P , l2 
2a i , m \vv\ ■ 
\u{q)\ 



1. Separable Ansatz 

The form of the gluon propagator A(p — q) appearing in Eq. ( |A7|) is not yet specified. 
Introducing the separable form in Eq. (|TTj ) and taking into account the symmetry properties 
of the functions /, W and U under p ■ P — » —p ■ P, which follow from charge conjugation 
symmetry, one obtains the following sets of integral equations for the scalar mesons: 

a. Scalar, C = + mesons 
f(p) = y / G(p 2 ) G(q 2 ) [T ff f(q) + T fW W(q) + T fU U(q)} , 

W(p) = ^ J -0^ p ■ q F(p 2 ) F(q 2 ) [T wf f(q) + T ww W(q) + T wu U(q)] , (A15) 
U(p) = | / -0^ p ■ q F(p 2 ) F(q 2 ) w(p, P) ■ u(<7, P) [H/,/ (?) + IW^(g) + !WC/(g)] 

6. Scalar, C = — mesons 

fb) = ffj0ji p ' 9 F{p2) F{q2) [Tfffiq) + T f wW{q)] ' 



W(p) = 3/^4 G(P 2 ) G(g 2 ) [T w/ /(g) + IWW(g)] , (A16) 
C/(p) = 0. 



The Tff, Tfw, ■ ■ ■ are given by Eqs. (|A8| , [A9|) . There are similar equations for the pseudoscalar 
mesons. 
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2. Form of the solution using the constrained separable Ansatz 



The separable form of the propagator causes the solutions of these equations to be 
proportional to the functions G and F. For the scalar, C = + meson, for instance, one finds 
that the solution is of the form 

f(p) = A/G(p 2 ), 
W{p) = -i\ w p- PF(p 2 ), (A17) 
U(P) = ^^(p 2 ) • 

Substituting the above Ansatz into Eqs.( [A15| ) yields a simple matrix equation of the 
form 





K(M) \ \ w , (A18) 



where K(M) is a 3 x 3 matrix whose elements are two-dimensional integrals that are com- 
pletely determined once ay and as are specified. This equation is then solved by adjusting 
the meson mass M until one of the eigen-values of K equals one. This procedure can be 
implemented by introducing an eigen- value, /i(M), on the left-hand-side of Eq.( |A18j ); solv- 
ing for /^(M) and the eigen-vector at each value of M; and repeating the process until one 
finds M such that fi(M) = 1. At this point one also has the Bethe-Salpeter amplitude for 
the bound state, which is characterised by the multiplet {A/, ^u}- 

Bound states of unequal mass quarks [see Eqs. (|l8l) and Sec. [HI A 2| , for example], are 



not characterised by a charge-conjugation quantum number, C. In this case the functions 
/, W and U are complex and Eq. ( |A17|) generalises to forms such as 



f(p) = \ lu G u (p 2 ) + \ ls G s (p 2 ) + q-P [\ 2 uF u (p 2 ) + \2 S F s (p 2 )\ +.... (A19) 
In this case the analogue of the matrix K(M) in Eq. (|A18|) is, in general, a 12 x 12 matrix, 



which reduces to a 10 x 10 matrix when residual symmetry under C is taken into account. 

APPENDIX B: BETHE-SALPETER AMPLITUDES IN SEPARABLE 

APPROXIMATION 

1. Scalar and Pseudoscalar Mesons. 

a. fi = u/d = f 2 



The most general form for the solutions of Eq. ( |33"D in the scalar and pseudoscalar 
channels are: 
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r^ d (p,p) = G u ( P 2 ) \ f i D -%\ wl .p 



-Xwp-Pj-P + Xu 



u{p)\ 



7 • u(p) 



T Vseu d ^ p)=Fu{p2 



iXfP ■P + \ w p-P>y-P-\ 



u 



u(p)\ 



7 ■ Hp) 



lib, 



i scalar 



(p, P) = i P u (p 2 ) X f p-PI D - G u (p 2 ) X wl -P 



(Bl) 
(B2) 

(B3) 
(B4) 



where P^ [P 2 = —1] is the direction-vector associated with P^, and u^{p,P) is defined in 
Eq. (|A4|). We note the covariants involving 7 • ■u(p) may be brought to a more familiar form 
through use of the identity 



1 



Hp) I 



7 • u{p) = 7 ■ p + p ■ P 7 ■ P. 



(B5) 



Solving Eq. (|3^) yields the following C = ±, pseudoscalar and scalar eigen- vectors: 

jpc -+ ++ — o+- 

X f 0.61 0.67 -0.52 0.11 
X w -0.045 -0.0075 0.084 0.26 
X v 0.0 -0.050 0.024 0.0 



(B6) 



which are normalised in accordance with Eq. (|T9|). 



b. /1 = u/d, f 2 



Solutions for the pseudoscalar Bethe-Salpeter amplitudes for the w-s-mesons are of the 
form 

r pseud (p,p) = 

\x lu G u (p 2 ) + X ls G s (p 2 )-p-P \X 2u F u (p 2 ) + X 2s F s (p 2 ) 

X 4u F u (p 2 ) + X 4s F s (p 2 

X 5u F u (p 2 ) + X 5s F s {p 2 ' 



(X3uG u (p 2 ) + X 3s G s (p 2 )-p-P 

1 

\u(p)\ 



i'j ■ P 



7 • w(p) > ^75 • 



(B7) 



The scalar amplitude has the same form but with ij 5 — > I D . The pseudoscalar amplitude 
given here corresponds to the K~ . The K + amplitude is obtained by making the replacement 
P» -> -Pit- 

The calculated pseudoscalar eigen-vector is 



J P X\ u Xu X 2u X 2s Xs u X 3s X^ u X4 S X^ u X§ s 
0- 263 390 -1.3 -3.3 -60 -97 2.9 7.0 2.8 6.7 



:bs) 



where each of the components is to be multiplied by 10 3 and the normalisation is in 
accordance with Eq. (|i~9"|) . No true scalar solution is found. 
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c. 7] Meson. 

The normalisation condition for the rjg p meson is 
d 4 k 



2P„ = N r . 



(27T)' 



■ x 



{|(cos^ P - ^sin^p) 2 [tr (T v (k, -P)d?S u (k + \P)Y v (k, P)S u (k - \P) 
+tr (T v (k, -P)S u (k + lP)T v (k, P)d?S u (k - \P)) 
cosOp + sm9 P ) 2 [tr (r v {k, -P)d?S s (k + \P)T v {k, P)S s {k - \P) 

+tr (r v (k, -P)s s (k + lP)r v (k, P)d?s s (k - |p))] } . 

The formula for the decay constant of the r]0 p meson is: 



(B9) 



P 2 U 



Nr. 



d A k 



■ x 



(2tt)4 

{|(cos^p - v^sintfpftr [7 ■ P^S u (p + \P)V n {p, P)S u (p - ±P) 
+ §(v/2 cosfl P + sin^p) 2 tr [7 • P l5 S s (p + \P)T v {p, P)S s (p - ±P)] } . 

The positive charge parity solution of Eq. (^) has the form: 

r„(g,P) = 



(BIO) 



\fuG u (q 2 ) + XfsG s (q 2 ) - X Wu G u {q 2 ) + X Ws G s (q 2 ) 17 • P 



lib- 



(Bll) 



Solving Eq. fl3"5|) leads to the following values of the the mass, decay constant and eigen- 
vector at the listed values of 6p: 



6 P 


5° 


0° 


-5° 


-10° 


-90° 


-95° 




0.549 


0.513 


0.475 


0.436 


0.357 


0.399 


fve P 


0.114 


0.111 


0.108 


0.105 


0.100 


0.102 


A/„ 


0.18 


0.23 


0.28 


0.33 


0.43 


0.38 


Xfs 


0.47 


0.41 


0.35 


0.29 


0.19 


0.25 


X\Vu 


-0.041 


-0.051 


-0.059 


-0.066 


-0.073 


-0.071 


-Ws 


-0.11 


-0.092 


-0.074 


-0.058 


-0.032 


-0.045 



(B12) 



The masses and decay constants are given in GeV. The last two columns correspond to the 
7/ flavour projection. 



2. Vector and Axial-vector Mesons. 

a. fi = u/d = f 2 

The on-shell constraint of Eq. ( |3"9"| ) entails that in constructing the general form of the 
vector meson Bethe-Salpeter amplitude we can work with the following transverse Euclidean 
covariants: 
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Pi, iL vll ■ P, Pll -Pp-P, IhtuvXpluPxPp 



(B13) 



The general amplitude is a linear combination of these covariants weighted by invariant 
amplitudes Ti{p 2 , P 2 , p ■ P). Employing the separable Ansatz, Eq. (|17"D, the Bethe-Salpeter 
equation, Eq. (|3"%|), cannot support contributions to T u (p,P) that are bilinear in p. Hence, 
at the mass-shell, the produced vector meson Bethe-Salpeter amplitude is 

I?(p, P) = plFuip 2 )^ + z^G u (p 2 )X 2 + ij 5 e„ vXpllJ ,p x P p F u (p 2 )X 3 . (B14) 

For the axial-vector meson the transverse Euclidean covariants are 

pllsP ■ P, 757J, Pvlsl ■ P, Pllsl -Pp-P, £nu\ P lnP\Pp ■ (B15) 

Again the terms bilinear in p do not contribute and, using the separable Ansatz, the produced 
axial-vector Bethe-Salpeter amplitude is 

Tl(p, P) = i^lGuip 2 )^ + ie pvXplp p x P p F u {p 2 )X 2 . (B16) 

As above, the Bethe-Salpeter equation is a matrix eigen-value problem. We obtain the 
solutions 

J PC Xi X 2 A3 

1— 0.075 -0.33 0.049 (B17) 
1++ 0.056 -0.28 0.0 

corresponding to the p/u and a^j f\ channels respectively. The Bethe-Salpeter amplitudes 
are normalised according to the vector and axial- vector generalisations of Eq. (|i~9D . 

b. fx = u/d, f 2 = s 

Consider the J p = 1~ K* + meson. As for the u meson, there are 5 transverse covariants, 
which can be taken to be those in Eq. (|B13|) except that the explicit factor of p ■ P is no 
longer necessary because w-s-states are not eigen-states of the charge conjugation operator, 
C. The general amplitude is a linear combination of these covariants weighted by invariant 
amplitudes ^{p 2 , P 2 ,p ■ P) where odd powers of p ■ P are allowed for the same reason. The 
separable Ansatz does not support contributions bilinear in p and hence, at the mass-shell, 
the produced K* + amplitude has the form 

r£(p, P) = p T v (f u (p 2 )Xiu + F s (p 2 )X ls ) + ill (G u (p 2 )X 2u + G s (p 2 )X 2s ) 

+ i-fip ■ P (F u (p 2 )X 3u + F s (p 2 )X 3s ) + ip T vl ■ P (F u {p 2 )X iu + F s {p 2 )X 4s ) 

+ i^e^xpluPxPp (Pu{p 2 )X 5u + F s {p 2 )X 5s ) . (B18) 

The amplitude for the i^*~-meson is obtained by reversing the sign of p, under which the 
kernel is invariant. The amplitudes for the J p = 1 + K\ meson states are simply 75 times 
the appropriate form of Eq. (pi8|) . 
The calculated eigen- vectors X\ are 
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J p 



Ai 
u 0.020 
s 0.046 
u 0.16 
s 0.34 



A 2 

-0.12 

-0.21 
-5.4 x 10~ 3 
4.0 x 10~ 3 



A 3 

-4.5 x 10~ 4 
-6.5 x 10" 4 
-6.7 x 10^ 3 
-1.5 x 10~ 2 



A 4 

2.0 x 10~ 4 
3.4 x 10~ 4 
-1.8 x 10~ 3 
9.9 x 10~ 3 



A 5 
0.014 
0.026 

6.8 x 10~ 4 

1.9 x 10~ 3 



(B19) 



c fi = s = f 2 

The Bethe-Salpeter amplitude for the J PC = 1 ss state [0] has the same form as 
Eq. QBTg) for the p/u but with G u -> G s and F u -> F a . That for the J PC = 1++ state 
[/i(1510)] is related in a similar way to Eq. (|B16|) for the a±/fi. The calculated eigen- 
vectors are 

J PC Ai A2 A3 

1— 0.049 -0.35 0.030 (B20) 
1++ 0.0044 -0.18 0.0 



3. Diquark Correlations 

All diquark eigen- vectors are normalised according to \\\ — 1- 



a. /1 = it/d = h 

The calculated eigen- vectors for the scalar (0 + ) and pseudoscalar (0~) diquark correla- 
tions are 

J P A 3 Afy Xf; 

0+ 0.96 -0.29 0.0 (B21) 
0- 0.15 0.58 -0.80 

The calculated eigen- vectors for the axial-vector (1 + ) and vector (1 _ ) diquark correlations 

are 

J p Ai A2 A3 

1- 0.12 -1.46 0.0 (B22) 
1+ 0.16 -0.98 0.11 



b. h = u/d, f 2 = s 
The calculated eigen- vector for the scalar (0 + ) diquark correlation is 

J P Ai„ Ai s \ 2u ^2s ^3u ^3s ^5s (Wl) 

0+ 498 802 -10.4 -33.2 -165 -282 5.2 14.9 4.7 11 { ' 
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where each component is to be multiplied by 10 3 . 

The calculated eigen-vectors for the 1 + and diquark correlations are 



J p 




Ai 


A 2 


A 3 


A 4 


A 5 




u 


168 


-4.7 


-0.029 


-7.0 


4.0 




s 


445 


6.3 


0.21 


-17.0 


39.6 


1+ 


u 


16.2 


-161 


-0.64 


0.13 


12.9 




s 


45.8 


-288 


-1.05 


0.21 


26.8 



where again each component is to be multiplied by 10 



c- fi = s = f 2 

The calculated eigen-vectors for the 1 + and V diquark correlations are 

J p Ai A2 A3 

I" 0.030 -1.64 0.0 (B25) 
1+ 0.090 -0.99 0.061 
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TABLES 



fi 


Jj 


%< 


a calc. Ge y2 


b c f f c GeV 2 


u/d 


0.482 


0.00811 


0.0413 


0.0281 


s 


0.580 


0.198 


0.0385 


0.0426 



TABLE I. Values of the fitting parameters, eg and rh^, used in constructing the constrained, 
separable Ansatz; and the values a/ ; and bf i , defined in Eqs. ( |i~5| ) and (|l^), calculated using them. 
The parameters eg and rh are chosen so as to fit = 137.5 MeV and f„ = 92.4 MeV; the 
parameters es and m s so as to fit mx = 493.6 MeV and fx = 113 MeV. The values of rhf % listed 
here correspond to m u u = 4.59 MeV and m s = 112 MeV. See Eq,([3~4|) and associated text for 
further details. 





mf} c - GeV 


_„calc. Dom 
m M 


Expt. 


7T (0-+) 


0.139 (fit) 


0.116 


7r ± (140),7r°(135) 


/o/ao (0 ++ ) 


0.715 


0.743 


/ (980)/a (982) 


o+- 


1.082 


1.092 


Not Seen 


o— " 


1.319 


1.299 


Not Seen 


K 


0.494 (fit) 


0.412 


#±(494), K° (498) 


K*o 


unbound 


unbound 


K *(1430) 


V(9p = 5°) 

v(0p = o°) 


0.549 
0.513 


0.472 
0.441 


7/(547) 


uj/p 


0.736 


0.755 


cj(782)/p(770) 


ai/fi 


1.34 


1.37 


01(1260)//! (1285) 


K* 


0.854 


0.866 


K*(892) 


K X 


1.39 


1.39 


#1(1270), #i(1400) 


<f> {ss l~) 


0.950 


0.957 


0(1020) 


ss 1+ 


1.60 


1.60 


/i(1510) 



TABLE II. Calculated meson masses compared with experimental values [25|, when known. 
The column labelled with the superscript "Dom" means that the quantity was calculated using only 
the leading Dirac amplitude; e.g., T n (p,P) cx i^G u {p 2 ) for the pseudoscalar; "unbound" means 
that in ladder approximation our constrained, separable Ansatz does not yield a stable bound state 
in the channel under consideration. 
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nf- GeV 


flf c - Dom GeV 


Expt. 


7T 


0.0924 (fit) 


0.056 


vr+(0.0924) 


K± 


0.113 (fit) 


0.76 


/if +(0.113) 


rj(9 = 5°) 
r,{9 = 0°) 


0.114 
0.111 


0.086 
0.082 


0.094±0.007 or 0.091±0.006 



TABLE III. Calculated weak decay constants compared with experimental values p5[ . The 
superscript "Dom" has the same meaning as in Table O. 



fl 


h 


J p 


M GeV 


M Dom Q e y 


u/d 


u/d 


0+ 


0.737 


0.653 


u/d 


u/d 


o - 


1.50 


1.52 


u/d 


s 


0+ 


0.882 


0.786 


u/d 


s 


o- 


unbound 


unbound 


u/d 


u/d 


1+ 


0.949 


0.958 


u/d 


u/d 


1- 


1.47 


1.48 


u/d 


s 


1+ 


1.05 


1.05 


u/d 


s 


1- 


1.53 


1.53 


s 


s 


1+ 


1.13 


1.13 


s 


s 


1- 


1.64 


1.64 



TABLE IV. Calculated diquark effective-masses. The superscript "Dom" has the same meaning 
as m Table 0. 
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